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Abstract 

We derive recursion relations for closed string correlation functions and scattering am- 
plitudes which hold to all orders in perturbation theory. 

1 Introduction 

Recent years have seen a renewal of interest in recursive calculation of the S-matrix. This 
has largely been in the context of gluon scattering using helicity variables [I] following 
[2] which is a more efficient and revealing calculational method than the usual LSZ and 
Feynman diagram expansion. At tree level these recurrence relations are completely 
understood, along with their connections to twistor theory [3]- [11] and there is now 
significant understanding of one loop results |12j-[22|. Recent discovery of an MHV 
action [23]- [25] opens the way for a systematic exploration of higher loop results. This 
program has lead to related results in QED, gravity and other theories [26]-|33] although 
results are again mostly limited to one loop. 

In a recent paper we presented recurrence relations between QED S-matrix elements 
[34] ( [35] for a similar approach to QCD) which displayed the same basic structure at 
tree level as the MHV recursion relations but extended simply and explicitly to all loop 
orders. In this paper we will show that the closed string S-matrix admits a similar 
recursive expansion using light cone string field theory [36]- [39]. The results will be 
valid to all orders in the loop expansion. 



2 Recursive expansion of the S-matrix 

The light cone closed string field action is 

where the string field <j) p+ is a function of light cone time r, the momentum p + conjugate 
to X- and the 24 transverse co-ordinates X. The measure is shorthand for 
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The operator h is the first quantised string Hamiltonian, 



We fix the length of the strings at 2n rather than the usual 2p + n. The latter convention 
naturally lends itself to explicit interaction calculations and though we always consider 
the interacting theory in this paper we will not need the precise details of the three 
string vertex. We therefore adopt the former convention which is more familiar from 
first quantisation. 

The generating functional of connected correlation functions, W[J, J*], is given by 
the functional integral 



o w[j,j^\/h 



= J 25(0,^) exp f^S[<j>,^] + jj DXdp+ 4 + <j> p+ - J p+ 4>1 
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This functional may be expanded in the number of sources J, and in powers of h, 
corresponding to a loop expansion, 

W[J, A = E / II D ( X <> Y i) 4 +1 ( X ^ • • • 4 +n (*n) J, + 1 (Xl) ■ ■ ■ J g+ JY m ) 

n,m— 1 i=l..n . . 

1=0 j=l. .m (4) 
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The subscript L indicates the loop order. Subsequent equations will become unwieldy 
unless we further condense our notation. When unambiguous we will suppress depen- 
dencies on the co-ordinates and momenta. We will also abbreviate 

( \T4> P+ (Xi) . . . P+ „ {X n )4 +1 (Fi) . . . 4>\ +m (Y m )\ > = (1 . . . n; 1 . . . m). (5) 

Then the expansion of the generating functional may be abbreviated to 

W\J,J*\ = Y [ (»")H) m Jt "jm h L {1 i m)L 
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Our goal is to evaluate the expectation value 



1 - ^t- 
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n 

which we may write in terms of functional derivatives with respect to the sources, 



We are not interested in free correlators so we separate W into 



W = J J f (V + i^-) J + W[J, J f ] 



(7) 



where W contains correlation functions of the interacting field. This decomposition 
reveals a divergence in the expectation value coming from the free part. This is because 
(f6j) involves operators evaluated at coincident spacetime points and therefore stands in 
need of regularisation. To do this we will insert the operator % into the free action 
acting on <^ , the effect of which is to shift the arguments of the field according to 

T £ 4 + [X J r]=</»t ++ JX + e 1 r + e]. 

where (e, e,e+) is some small change in position and momentum. This removes the 
divergent term. Evaluating (|6]) and letting e — + we find 
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The effect of the final two terms is simply to count the number of string fields so that 
the 4>"j 4> im correlator in the perturbative expansion of W is multiplied by — (n + rn). 
Including the regularisation we may also write the expectation value as 



C=0 



where the action S is defined by 



h 



s = j<& + [id T - — J (i + CT e )0 p+ + s 3 [4>t + , 4> p+ ]. 

The functional integral in ([9]) is given by the usual Feynman diagram expansion of the 
generating functional except the propagator G is replaced by 

The inverse operator may be defined as a formal power series in f. The derivative, as £ 
goes to zero, multiplies the exponential of W/h by a sum of Feynman diagrams derived 
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from the usual sum over connected diagrams as follows. Each connected diagram is 
separated into E + I similar diagrams, where E (I) is the number of external (internal) 
lines. In each of these diagrams the £ derivative acts on one propagator. As £ goes to 
zero this propagator, G, is replaced by —T,_G. 

As we let the regulator e — + we therefore recover the sum over connected diagrams 
where each is multiplied by a factor of — E — I. Since the interaction is cubic the number 
of internal lines in a connected diagram is a function of the number of external lines 
and the loop order L given by /f := E + 3(L — 1). The expectation value may then be 
written 



y r (t")( i) m Jt n jmh l {1 1 _ m) (n m r + m) 

t—1 ] ram 
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Equating the expressions (J5J and (|10[) we find 



/ | , J J ft (l...n;l...m) L 7 ^ 

n,m,L 



L 



5W/ Q . ft \<W t /• 5 . ft \<W 



(11) 



The free part of the generating functional has dropped out of our equations. Equating 
by order in J, and ft we find 

l2 +m {l...n;l...m) L = J T)X (d-r + 1^-^(1 .. .n, X;l .. .m,Y) L _ : 
+ /^^C K/r<l...n-r;l... S ,X) I/ (V+i 



MX) 

2p+ 



(1 . . . r, X; 1 . . . m — s) 



L — L' 



r=0..n 
s — 0. .m 
L'=0..L 



(12) 



where i?"' s m = l/(r!s!(n — r!)(m — s!). The sum over a is over all distributions of 
the external data {X r ,Y r } between the two correlation functions. We claim that these 
relations are recursive in the number of string field operators and in the loop order, but 
the initial and final terms of the sums may appear to spoil this. On the left hand side 
is a correlator of n + m fields at loop level L, while the right hand side is built from 
pairs of correlators nearly all of which are for < n fields at < L loops. However when 
r = and s = m — 1 one correlator contains a total of m + n fields and may be of 
loop level L. Observe though that this is coupled to a correlator of two fields at loop 
level 0, which is just the free propagator. By construction this object does not appear 
in W and this term vanishes from the sum. Therefore, although the right hand side 
contains correlators of n + m fields they are of loop order at most L — 1. When r = n 
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and s = the right hand side contains one field correlators (tadpoles and higher loop 
contributions) coupled ton + m + 1 field correlators. When the latter is of loop order L 
the one point function is tree level. Since there are no tree level tadpoles, by definition, 
this term vanishes and correlators of n + m + 1 fields are of loop order at most L — 1. 
Equation (I12|l therefore gives a genuinely recursive expansion of the partition function 
in the number of correlated fields and the loop order. 

It is now a short step to go to a similar expansion of the S-matrix. One method is 
to repeat the above analysis with the sources J and J replaced by J'(d T + ih/2p+) 
and J(—d T + ih/2p+). Differentiation with respect to the Fourier modes of J and 
brings down the Fourier transform of correlation functions with amputated external 
legs, i.e. (generally off-shell) S-matrix elements. The calculation proceeds as before. 
Alternatively, we may directly amputate and Fourier transform equation (|12[) . 

The creation and annihilation modes may be extracted from the string field in the 
usual way (see appendix) , 

4i n p,0} = / DX /{0(X)e- lE - T+lp ' x ^ + P0, 

^+,p,o} =/dx / {0 (XK E -— p - x P+ (X). 

The LSZ reduction formula for S-matrix elements follows using the usual methods. De- 
noting external data {p+ r , E-, p r } by simply r, the S-matrix element 

5(E- in - £-out)5(Pin - Pout)<5(jJ + in ~ P+out)A(l . . . n\l . . . m) L 

is given by a product of terms 



for each 'out' state and 



(14) 



J dr r J DX r e* B - (l r>Vr>P+r)-r-Pr~r) f{i} ^ Xr) U T + 

' state and 

J dr; I BY r e -^-^, q+r )<-^^) f{i}Mr) ^ dr , + iMXlT) 

for each 'in' state, all acting on the correlator 

(o|t^ + (Xx) . . . ^, +n (Jr B )^ +1 (yo . ..4>t + JY m )\o) L . 

We apply these operators to (|12p . For the integrated (internal line) data appearing in 
the correlation functions we may use the spectral decomposition of the propagator in 
order to introduce the prefactors (|13[) and ()14|) associated to the off-shell string fields; 

G fXrYr'l-tf^-^lT /" d 24 qdg_ (T _ r > )+iq(x _ y) / W (X)/ W (Y) 

(15) 
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We then find that the scattering of m to n strings is given by 



r L +m A(l...n\l..,m) L =J2 




id^k 



2k + 



A(l...n, {k, {l}}\l...m,{k,{l}}) L 



-i 



n-r\l...s,{q y {l}}) L 



l 



_ h(q,{Q) 

2q + 



x A{1 . . . r, {q, {1}}\1 . . . m - s) L _ L , . 



(16) 



In the first line the measure is cPfc = (27r) _25 d 24 kdfc+ dfc_. In the second line the sum 
is as in (|12[) but now includes a sum over all possible excitations {/} of the string. The 
momentum q — {g+,q, g_} is forced, by the overall delta functions in the scattering 
amplitudes, to conserve momentum between pairs of amplitudes. If we restrict to tree 
level then the integrated terms in (|12p and (|16|l go out. Beyond tree level these terms 
remain and describe how L-loop amplitudes of E external strings receive contributions 
from L — 1 loop amplitudes of E + 2 external strings - the two additional external legs 
are sewn together using a propagator to form a closed loop. 



We have derived recurrence relations between scattering amplitudes of closed strings 
using light cone string field theory. We emphasise that in contrast to other approaches 
our recursion relations hold explicitly to all loop orders. 

In one sense this recursive structure is more natural in string theories than in particle 
theories where worldsheets may be joined together to form Feynman diagrams using Car- 
lip's sewing method [40 • Two worldsheets are sewn together by acting on two boundaries 
with the first quantised Hamiltonian (which gives the correct moduli space measure on 
the sewn worldsheet), inserting a closed string propagator between the boundaries and 
integrating over the common boundary data. Our results are effectively an application 
of this method, albeit in a Lagrangian formalism; pairs of amplitudes are sewn together 
using the inverse of the two point function, which is the inverse of the first quantised 
Hamiltonian. 

We have not needed the explicit form of the three string interaction. Any interaction 
of the form 



with V some integral kernel (such as the usual delta functional in all variables) will lead 
to the same recursion relations - we may in fact any number of cubic terms with any 
number of derivatives acting on the fields, provided we satisfy physical constraints such 
as reparametrisation invariance. The chosen interaction kernel provides the initial data 



3 Discussion 




G 



for constructing amplitudes recursively. We note that our method extends to non-flat 
spacetimes, for example the Plane Wave background ■ The recurrence relations hold 
as in (I12[l and (|16p . with h replaced by the plane wave Hamiltonian, 

hp W = ^flat + J da ^i 2 p+X(cr) 2 . 

o 

Formally, our arguments appear to extend immediately to any supersymmetric string 
field theory with a cubic vertex, e.g. the Heterotic and Type IIAB strings - the dif- 
ferences would amount to replacing the quadratic form in (|12p and (|16p with that ap- 
pearing in the superstring action. This would remove the tachyon divergences of the 
bosonic string which we have not treated here. Superstring field theory suffers however 
from contact divergences [42] [43] [44] which are absent in the bosonic case. The reg- 
ularisation of these divergences generally requires the addition of counter terms which 
are of higher than cubic order. These higher order terms would considerably complicate 
the form of the recursion relations we have derived. An interesting question for future 
study is whether is it possible to regulate contact term divergences while maintaining 
the simple structure of the bosonic relations. 



A Conventions 

For completeness we present our conventions. The equations of motion which follow 
from the action |l} are 

{idr - 7T— h)<P P+ = 0, (-id T - i^-h)4>l + = 0. 
To solve these equations we first expand the transverse co-ordinates in a suitable basis, 

^ oo 

X(cr) = — == H ■= > x n cos na + x n sin no, 



^ oo 

7^ ^ ' 



5 1 d 1 ^ q q 

■ cos na + sin na. 



v n — 1 



This allows us to expand h in modes such that the equations of motion become partial 
differential equations which we solve by separation of variables. With a slight abuse of 
notation, the separated equations are 



If d 2 



id T (f>p + 








dl<l>P+ 


2 J 

= -P <Pp+, 










for % = 1 . 


..24, 


) 2 W 




for % = 1 . 


..24, 



(17) 
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with the constraint 

E^ = ^(~p 2 -Y t E n ,i + E n \ (18) 

We recognise (|17|l as Hermite's equation, for well behaved solutions to which we must 
have 2E n /n = 21 + 1 for I a non- negative integer which we call the level. The solution 
of level I is 

H lnt {V^^)e' n{ < )2 /2 , (19) 

where Hi is a Hermite polynomial. This means that for each mode of the string, de- 
scribing a quantum harmonic oscillator of frequency n, in each dimension transverse i, 
there are an infinite number of solutions to the equations of motion representing the 
infinite number of excitations of the oscillator, labelled by The constraint equation 
(|18l) becomes 

24 oo 

2p+E- -p 2 - ^^nl n ,i +nl n ,i + n = 0. (20) 

i — 1 n — 1 

The sum over the final term diverges, and must be regulated. Zeta function regularisation 
implies 

E- = ^-(-2 + Y j nln.i + nTn.il = ~ + o ~ ~ 2 > ( 21 ) 

which is proportional to the first quantised string Hamiltonian Lo + -Lo which treats Z as 
the number of excited string oscillators of frequency n. The full solution to the equations 
of motion is then 

„ P+ ,r,x) = y - 

where the functions /{;} are defined by 

The product over n and i is the product over the separated Hermite solutions for each 
string mode x n in dimension i. There is then a sum over all possible combinations 
of assigning level values to each of these solutions, representing the infinite number of 
independent solutions to the equations of motion. We have included a normalisation 
factor in the Hermite polynomials, which obey 

J Ax e~ nx2 H h (V7ix)H l2 (^x) = 2 l l\{n /n) l/2 8 hM . (24) 
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